Random walks and diameter of finite scale-free networks 
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Dynamical scalings for the end-to-end distance R ee and the number of distinct visited nodes N v of 
random walks (RWs) on finite scale-free networks (SFNs) are studied numerically. (R ee ) shows the 
dynamical scaling behavior (Ree(i,t)) = I ( r y,N)g(t/£ ), where £ is the average minimum distance 
between all possible pairs of nodes in the network, N is the number of nodes, 7 is the degree exponent 
of the SFN and t is the step number of RWs. Especially, (R ee (I,t)} in the limit t — * 00 satisfies 

the relation (R ee ) ~ f ~ d a , where d is the diameter of network with d(i) ~ In N for 7 > 3 or 
d(£) ~ In In TV for 7 < 3. Based on the scaling relation (R ee ), we also find that the scaling behavior 
of the diameter of networks can be measured very efficiently by using RWs. 

PACS numbers: 89.75.Hc,05.40.Fb,89.75.Fb 



For almost a decade there have been many studies 
on the topological properties of complex networks, since 
many structures of physically interacting systems are 
shown to form nontrivial complex structures [H, [2] . In 
these studies much effort has been put to investigate the 
physical origin of complex networks 0, H[ . It has been 
found that most of real web-like systems share several 
prominent structural features: scale-free degree distri- 
bution, high clustering, and small average path length 
(APL) P, E II II- It is well known @, [f that random 
networks (RNs) [8[ and small-world networks (SWNs) 
9} have very small APL or diameter d which scales as 
d ~ In A, where A is the number of nodes. Recently, 
Cohen and Havlin [l(| studied the diameter of scale-free 
networks (SFNs) Q whose degree distribution P(k) satis- 
fies a power law P(k) ~ A: -7 . Using analytical arguments 
they showed that 



d(N) 



In In A 

In A/ In In A 

In A 
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7 = 3 
7 > 3. 



(1) 



For measuring the APL, the breath-first algorithm (BFA) 
[Tl| are mainly used. This algorithm is known to scale 
as 0(N 2 ). In this paper, we will show that random walk 
(RW) can provide more efficient method to find the scal- 
ing behavior of d like Eq. (QJ. 

The early studies on complex networks mainly focused 
on their topological properties. Recently, the physical 
systems whose elements interact along the links in com- 
plex networks have drawn much attention. Furthermore 
a number of studies have focused on the effects of the un- 
derlying topologies on the dynamical properties of such 
systems. Many dynamical systems on complex networks 
show rich behaviors which are far from the mean-field 
expectations and affected by the underlying topology 
For example, the dynamical properties of RWs 
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on complex networks have been shown to be cl osely re- 
lated to the topology of underlying networks fl4l.lTa[l6l|. 
Especially, the average number of distinct visited sites 
and the average end-to-end distance of RWs on SWNs 
are known to satisfy the scaling law fl6( ] 



Q(p,N,t) =O sat F(p%pN) 



(2) 



Here, p is a density of shortcut and O sat is a saturated 
quantity. 

It is well known that the SFNs have heterogeneous 
structures in which nodes with anomalously large number 
of degree and nodes with small degree co-exist [3, 3 . In 
SFNs, the dynamical properties of several systems are af- 
fected by the second moment of degree distribution (A: 2 ) 
[13, EH [S H3| ■ It is interesting to study how the struc- 
tural heterogeneity affects the scaling properties of RWs 
as changing the degree exponent 7. In this paper, we 
mainly investigate the dynamical scaling relation for the 
end-to-end distance R ee of RWs on finite SFNs with var- 
ious 7. From the scaling relation, it will be shown that 
RWs on SFNs can provide much more efficient method 
to measure the scaling behavior of the diameter of finite 
SFNs than BFA. For the complementary purpose we also 
study the scaling relation for the number of distinct vis- 
ited nodes N v . 

All the models are defined on a graph Gr = {A, AT}, 
where A is the number of nodes in the network and K 
is the total number of degrees with the average degree 
(k) = 2K/N. Networks of the sizes A = 10 3 ~ 10 6 are 
used. To generate SFNs, we use the static network model 
introduced by Ref. [2l[ with (k) = 4. Initially a random 
walker is placed on a randomly chosen node denoted by s. 
At the time step t + 1, the walker jumps to a randomly 
chosen node among the nearest neighbor nodes of the 
node where the walker is at t. The probability P(i, t; s, 0) 
to find the walker at node i at t thus follows the relation 



P(i,t + l; S ,0) = ^4^ P (^ s <°) 



(3) 
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with the condition P(i, 0; s, 0) = 5i, s [1 
adjacency matrix whose elements are Aij = 1(0) if two 



. Here, Aij is the 



2 



nodes i and j are connected (disconnected). All quanti- 
ties are averaged over 100 network realizations and 1000 
different initial positions of RWs for each network real- 
ization. 

First, let us discuss the average end-to-end distance 
(Ree) of RW on SFNs. At each time step i, we measure 
the shortest distance, R ee (t), from the node where the 
random walker is to the node s. By averaging i? ee (i) 
over different initial positions and network realizations 
we get (R ee )- Here the shortest distance between two 
nodes in networks means the shortest path length or the 
minimal number of steps between them. 

On a D-dimcnsional regular lattice, the average end- 
to-end distance follows (R ee ) ~ t v with i/ = 1/2 [22I ]. 
However, on finite systems, (R ee ) saturates to a finite 
value in the limit t — > oo. In this limit, the probability 
for each site to be occupied by a random walker is the 
same for all the sites and thus (R ee ) ~ I, where £ is 
the minimum distance averaged over all possible pairs 
of sites [16]]. On a SWN, for small t (<C 1/p 2 ) regime 
(-Ree) again satisfies the same relation as on the regular 
lattice ((R ee ) ~ t 1 / 2 ), because the random walker cannot 
actually see shortcuts. But in the large t regime (t ^> 

l 2 /p), (R ee ) ~ £ ■ For the intermediate time region 



1/p 2 <C t -C v £ 2 /p, Almaas et al. expected that (R ee ) ~ 
£ satisfies the relation (R ee ) ~ tp jig ]. 

A typical dependence of (R ee ) on t for finite SFNs is 
shown in Fig. [2 Figure 1 shows the measured (R ee ) 
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FIG. 1: Average end-to-end distance, (R ee ) on SFNs with 
7 = 2.4 for TV = 10 7 . (a) Plot of (Ree) against t in log- 
log scale. The solid line correspond to u = 0.90. The inset 
displays the dependence of the measured v on 1/N. From the 
inset, we expect that v — * 1 in the limit TV — * oo. (b) Plot of 
(Ree{£,t — » oo)) against £ in log-log scale with £ ~ In TV. The 
solid line corresponds to a = 0.93. 

on SFNs with 7 = 4.3 for N = 10 7 . From the early-t 
behavior of (R ee ) (or the data for t < 10 in Fig. [2(a)) 
the obtained value of the exponent v for the relation 

(R ee (t))~t v (4) 
is v ~ 0.90(5). The inset in Fig. QJa) displays the depen- 



dence of the measured v on the network size N, which in- 
dicates that the values of v approach to 1 as increasing N. 
This result agrees with the expectation on the SWN (IH ] 
when p — > 1. The (R ee (t)) does not increase indefinitely, 
but rapidly reaches a saturation value (R ee (£,t — > 00)^ 
after a very short crossover time r ee . Since r ee is nearly 
equal to 10 or slightly larger than 10 even for the very 
large network size (or N = 10 7 ), the expectation of v = 1 
from the data for nearly one decade of t or so seems to 
be physically dangerous. However, since an earlier ana- 
lytical study on the walks on a Caylcy tree also suggests 
such v = 1 behavior [H, v — ► 1 as N — > 00 is physically 
plausible one. We have also measured v for the size upto 
N = 10 7 and obtained the similar data as in Fig. [lja). 
We thus expect that v — > 1 as N — ► 00 regardless of 7. 

In Fig. HKb) we display the (R e e(l t -> 00)) for 7 = 4.3 
as a function of £ with £ ~ In N. As shown in Fig. [2(b) 
(Ree{£, t — > 00)) for 7 = 4.3 satisfies the power-law: 



(R ee (£,t ^ oc)) ~£ a . 



(5) 



with a = 0.93(2). <^i? ee (^, f — > 00)) for various 7 is also 
found to satisfy the power law ([S]) very well. The ob- 
tained a's for various 7 by assuming £ ~ In In N for 7 < 3 
and £ ~ In N for 7 > 3 are displayed in Table Q] 

From the Eqs. (01 and the dynamical scaling re- 
lation 



(Ree(£,t))^T( 1 ,N)g(t/f). 



(6) 



is expected. The scaling function then g(x) satisfies the 
relation 



9(x) 



X v , X < 1 

const. , x ^S> 1. 



(7) 



The dynamical scaling relation ^ also physically means 
that T ee scales as r ee ~ £ z ■ The dynamic exponent 2 can 
thus be evaluated from the measured a and z through the 
relation z — a/v (see Table Q|. Fig. [5] shows the scaling 
plot of (-Ree) measured on SFNs with 7 = 2.4, 3.0, 4.3 
and on random networks (RNs) for N = 10 3 , 10 4 , 10 5 and 
10 6 using Eq. ([6]). As shown in Fig. [2l (i? ee )'s for various 
N collapse very well into a single scaling curve with the 
exponents listed in Table [I] for each network topology. 
Note that £ in networks represents the average diameter 
by definition [l[. From the scaling collapse, we find the 
dependency of I on 7 and N: 



£(j,N) 



In In TV 
In A 



2 < 7 < 3 
7 > 3. 



(8) 



For SFN with 7 = 3, it is known that the diameter 
scales as d ~ In N/ (In In N) in the Barabasi- Albert (B A) 
model [lil with looped structures [1, [25[ . In contrast the 
BA network without looped structures or tree-type BA 
network satisfies d ~ In N 0, [25[ . For the best scaling 
collapse we use the relation d ~ In A for 7 = 3 even 
though the static SFN [H| has loops (see Fig. ^b)). To 
check the scaling behavior d ~ In N for 7 = 3 obtained 
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FIG. 2: Scaling plot of (R ee ) on SFNs with (a) 7 = 2.4 , (b) 
3.0, (c) 4.3 and on (d) RNs for AT = 10 3 , 10 4 , 10 5 and 10 6 . 
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FIG. 3: Dependence of the diameter of static SFN with 7 = 3 
on N. (a) The relation d ~ IniV is well satisfied, (b) The 
relation d ~ In N / In In JV deviates when N increases. 



by RW method, the diameter of SFN for 7 = 3 is directly 
measured by the BFA method (see Fig. [3]). As can be 
seen from Fig. (3]^a), the relation d ~ In N is well satisfied 
for the static SFN with 7 = 3 and agrees with the result 
obtained by RW method or from Fig. [5J However, the 
relation d ~ In N / In In N deviates when N increases (Fig. 

Mb))- 

The results from Eqs. (6) and (7) provide another 
way to find the scaling behavior of d in SFN. The result 
implies that the computing time needed for the measure- 
ment of scaling behavior of d by RW method increases as 
0((ln N) z ) or 0((ln In N) z ). Since the computing time to 
measure d by BFA increases as 0(N 2 ), the RW method 
is far more efficient to find scaling behavior of d. 

For comparison's sake we now explain the average 
number of distinct visited sites N v (t) on SFNs with var- 



TABLE I: Estimates of the exponents a, z for SFNs with 7 
and random network (RN). (z — a/u) 
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FIG. 4: This figure shows the scaling collapse for N v on SFN 
with 7 = 4.3 and 2.4 (inset). 



ious 7. Fig. [4]shows the scaling plot of N v (t) for 7 = 4.3 
and 2.4 against t/N for N = 10 3 , 10 4 , 10 5 and 10 6 . N v (t) 
on SFN satisfies the scaling relation 



N v - Nf(t/N) 



with 



1 , 1 « 1 

1 , X > 1. 



(9) 



(10) 



We have checked that the scaling relation ^ holds very 
well for SFNs with any 7. 

The scaling relation ^ is slightly different from that 
of the regular lattice and the SWN. On the infinite D- 
dimensional lattices N v depends on D. In the limit of 
t -> 00 [II], N v - Vt in D = 1, N v - t/\nt in D = 2 
and N v ~ t for D > 2. On SWNs, Almaas et el. 
showed that N v ~ y/t for t < ^ 2 (with £ = l/p),N v ~t 
for ^ 2 < t < N£_ and N v = N for t > N£. This means 
that if the random walker does not reach the shortcut, 
then the walker always sees the regular lattice structure 
(D = 1). However, if the walker meets the shortcut, then 
the behavior of walker follows the mean-field result and 
finally N v saturates to N due to the finite-size of the 
network. However, each node in SFN can be connected 
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to any other nodes and SFN can be regarded as an in- 
finite dimensional structure. Thus, N v on SFN follows 
the mean-field behavior after the first few steps. From 
Eqs. © and (jTUJ) the time t x at which N v = N is given 
by t x — N. This implies that the walker can sample 
a new region at each time step by using shortcuts until 
it visits all nodes in the network. Thus the statistical 
properties such as the average of certain quantity can be 
significantly enhanced even for small RW steps. 

Moreover, since (R ee (t)) ~ i 1/2 and (N v (t)) ~ t for 
D > 2, we can obtain a relation between the N v and the 
number of accessible nodes (or volume) N ac within the 
radius (R ee (t)) for D > 2: 

Thus, N ac /N v diverges as t — > oo which indicates the 
transient behavior of RWs [27| . The transient behavior 
becomes much more striking in the limit D — > oo. From 
r ee ~ (In TV) 2 ( or (In In TV) 2 and t x ~ N, we know r ee < 



t y on SFN. This result means that the walker effectively 
moves to the end of SFN without visiting all nodes. As 
a result, we can measure the scaling behavior of d(N) in 
much shorter time than BFA algorithm (for example, see 
Figs. Hand [3}. 

We investigate the scaling properties of RWs on SFNs. 
We measure the end-to-end distance R ee and the number 
of distinct visited sites N v . From the scaling ansatz for 
(R e e) as Eq. ([6]), we have measured the scaling exponents 
a, v and z for (R ee ) on various networks. From the 
scaling relation we find the dependency of £ on 7 and TV. 
Based on the simple scaling arguments, we also discuss 
the theoretical reasons why scaling behavior of d(N) of 
networks can be measured far more efficiently by RW 
method than by BFA. 
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